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A proposed shear sensor calibrator consists of a rotating wheel with the sensor mounted 
tangential to the rim and positioned in close proximity to the rim. The shear stress 
generated by the flow at the sensor position is simply r, v = prw/h, where pis the viscosity of 
the ambient gas, r the wheel radius, co the angular velocity of the wheel, and It the width of 
the gap between the wheel rim and the sensor. With numerical values of p= 31 pPa s (neon 
at room temperature), r = 0.5 m, co = 754 s 1 (7200 rpm), and h = 50.8 pm, a shear stress of 
t„= 231 Pa can be generated. An analysis based on one-dimensional flow, with the flow 
velocity having only an angular component as a function of the axial and radial coordinates, 
yields corrections to the above simple formula for the curvature of the wheel, flatness of the 
sensor, and finite width of the wheel. It is assumed that the sensor mount contains a trough 
(sidewalls) to render a velocity release boundary condition at the edges of the rim. The 
Taylor number under maximum flow conditions is found to be 62.3, sufficiently low to 
obviate flow instability. The fact that the parameters entering into the evaluation of the 
shear stress can be measured to high accuracy with well-defined uncertainties makes the 
proposed calibrator suitable for a physical standard for shear stress calibration. 
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Nomenclature 

radius of sensor 

Fourier coefficient 

area of sector of wheel 

correction for radius of wheel 

correction for finite width of wheel 

correction for flatness of sensor 

g(x) gap at location x 

gap at location x = 0 

gap in region beyond sensor 

wave number 

modified Bessel Function, first kind 

modified Bessel Function, second kind 

Mach number 

radius of wheel 

radius of trough 

Taylor number 

flow velocity 

r-component of flow velocity 
6 -component of flow velocity 
z-component of flow velocity 
peripheral speed of wheel 
width of wheel 
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coordinate in plane of sensor normal to axis of wheel 

coordinate along axis of wheel 

angular coordinate 

absolute viscosity 

density 

shear stress for rectangular coordinate system 

shear stress in cylindrical coordinate system 

shear stress based on a minimum rotation rate of 60 rpm 

shear stress at a maximum flow rate limited to M= 0.29 

shear stress at a maximum flow rate of M= 0.84 

parallel-plate shear stress 

angular velocity 
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An overbar indicates a mean value. 


I. Introduction 

I N order to calibrate or characterize a shear stress sensor, the sensor must be exposed to a known shear stress. With 
the advent of micro-electrical mechanical systems (MEMS) shear stress sensors the size of new sensors is such 
that new calibration techniques become feasible. There are two methods commonly used to calibrate these sensors. 
The first is to infer the shear stress in-situ using some other measurement technique. This methodology typically 
involves the use of Preston or Stanton tubes and can only be done under certain limited conditions but is useful for 
fully developed turbulent boundary layer research. The other method used is to create a known flow thereby 
generating a known shear stress. To generate a known shear stress a precisely defined flow must be generated. 

A. Existing Methods for Generating a Known Shear Stress 

1. Pressure Driven Flows 

One example of a pressure driven flow is a flow channel. A flow channel or laminar flow cell consists of two flat 
parallel plates of a finite width with the test gas, typically air or nitrogen, flowing between the plates and is one of 
the most commonly used methods for static shear stress calibrations. As long as the flow is laminar the velocity 
profile is well known and can be written in the following form: 


h 


(1) 


where /( y) is the velocity in the direction of the flow, y is the normal distance from the centerline of the flow 
channel, U-, is the centerline velocity, and h is one-half the channel height. From this equation the shear stress in this 
type of channel is found to be (ignoring compressibility effects) and assuming adiabatic wall conditions: 


r 


w 



(2) 


where the negative sign is an artifact of the coordinate system. Additional corrections for compressibility have also 
been developed. This analysis only applies for static shear stress measurements. One disadvantage of this flow is 
that the velocity profile is parabolic while the laminar sub-layer has a linear velocity profile. This may not be an 
issue for sensors that directly detect the shear stress (e.g., floating element sensors) but any method that infers the 
shear stress (e.g., thermal sensors) must determine if there is sensitivity to the shape of the velocity profile. 

2. Displacement Driven Flows 

This category refers to plane-wave tubes with either a circular or rectangular cross section. These devices rely on 
Stokes layer excitation to generate a boundary layer that is dependent on the excitation frequency. This device is 
used for low level dynamic characterization. Versions of this concept have been reported in the literature for both 
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zero mean flow and non-zero mean flow. Analysis of these devices is more complex than that of the simple analysis 
of a flow channel shown earlier and may be found in Ref. 1 and Ref. 2. 

3. Boundary Driven Flows 

The proposed device is a new design for a boundary driven flow. A boundary driven flow is one where one or 
more of the boundaries moves in a continuous manner and the flow is generated by the no slip condition at the 
interface between the moving surface and the test media. One such device uses the face of a rotating wheel to 
generate a known flow with a linear velocity profile. This flow has components in the r, 0, and z direction. This is a 
complicated flow field even though it has a linear velocity profile. 3,4 

B. Motivation for Development of a Rotary Flow Channel 

The primary motivation for designing a rotary flow channel is to enable precision calibration of MEMS shear 
stress sensors. Precision calibration requires a method for generating a known shear stress where the shear stress is 
dependent on a limited number of measurable or controllable quantities that can be determined with a known 
accuracy. The rotary flow channel has the potential to meet the requirements for a primary calibration standard for 
shear stress sensors. This device is intended to be one aspect of a comprehensive calibration scheme for MEMS 
shear stress sensors. 


II. Principle of Operation 

Figure 1 shows two parallel plates, the first moving at a constant speed U with respect to the second. The 
velocity profile is linear between the plates. The shear stress r 0) on the surface of the stationary plate is found to be 


du U 

r w =n— = n — , 
Sy g 


( 3 ) 


where // is the absolute viscosity of the medium, u the x-componcnt of velocity, x the direction of motion, y the 
normal direction, and g the gap between the plates. Equation (1) suggests, in principle, an accurate method of 
generating a known shear stress for the calibration of skin-friction sensors, since the quantities on the right hand side 
can be determined very accurately. However, a calibrator based on rectilinear motion of continuous duration will 
prove impractical. (Even a 
device based on a tape- 
drive mechanism has 
severe practical limitations, 
especially regarding the 
stability of the tape 
motion). If the tape 
fabrication issues can be 
resolved we would 
recommend that a device 
based on this design be 
fully analyzed and 
considered for use as a 
device for the calibration of 
shear stress sensors. 

However, these considerations suggest that a more practical approach is to utilize a rotary motion for the generation 
of the flow needed to produce the shear stress. 

Such a rotary flow channel is illustrated in Fig. 2. It consists of a wheel, sensor mount, trough, test medium, and 
supporting instrumentation. The wheel, having a radius R\, rotates at a fixed angular velocity a>, corresponding to a 
speed U=R\OJ at the surface of the wheel. The sensor mount, a flat plate containing the flush-mounted sensor to be 
calibrated, is positioned to realize a prescribed gap g. The trough has three functions: to entrain the ambient medium 
for flow into the gap, to suppress an axial component of the velocity, and to insure a velocity-release boundary 
condition at the edges of the wheel. Multiple configurations of the side wall design of the flow channel were 
considered and this version was selected for both ease of manufacturing as well as analytical simplicity. The 
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Figure 1. Velocity profile between parallel plates, top plate 
moving and bottom plate stationary. 



supporting instrumentation monitors the wheel speed, gap, and temperature and pressure in the gap. The use of a 



r 

a) b) 

Figure 2. Velocity profile between parallel plates, top plate moving and bottom plate stationary. 
Sidewalls of trough omitted from a) for clarity. 

gas other than air, as well as safety considerations, dictates that the entire calibration assembly be situated in a 
robust, sealed test chamber. 

Although the test medium can be air, the use of neon has several advantages. First, since it has the highest 
absolute viscosity of any common gas, a given shear stress can be generated at a lower wheel speed and larger gap 
than for air or other gases. Secondly, the high sound speed permits the generation of a given shear stress at a lower 
Mach number - an important consideration with regard to stability and compressibility effects. Thirdly, the high 
thermal conductivity facilitates the cooling of the medium and helps to maintain isothermal conditions in the gap. A 
comparison of pertinent physical properties between air and neon is shown in Table 1. 

Table 1: Properties of Air and Neon at 20°C. 



Air 

Neon 

Absolute viscosity p, Pa*s 

1.83E-05 

3.11E-05 

Speed of sound, m/s 

343 

449 

Thermal conductivity, W/(m*K) 

56.2 

107 

Density p at 1 atmosphere, kg/m 3 

1.21 

0.839 
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The levels of attainable shear 
stress can be substantial. If we 
ignore the effects of compressibility 
and the possibility for violating the 
no slip boundary condition, consider 
a wheel having a radius R 1 = 0.5 m, 
rotating at m = 7200 rpm 

(corresponding to U = 377 m/s) in 
neon, with a gap g x = 76 pm. The 
shear stress generated at the surface 
of the trough is then 154 Pa. 

III. Analysis 

A. Shear Stress Without 
Corrections 

For a rotary flow channel, 
corrections to Eq. (3) must be made 
for the finite radius and width of the 
wheel and for the flatness of an 
inserted sensor. The first part of the 
analysis will consider a trough of 
circular contour; then a correction 
will be made for the flatness of the 
sensor. 

Let the cylindrical coordinate 
system ( r , 6, z) have its origin at the 
center of the wheel, as is shown in 
Fig. 3. We assume that the velocity 
components ug = ug (r,z), u r = u~=0. 
The boundary conditions are: 



-w w 
2 1 2 

Figure 3. Cylindrical coordinate system for 
analysis of rotary flow channel showing both front 
and side views. 


u g (R„z) 

= U 

(4a) 

u g (R 2 ,z) 

1 = 0 

(4b) 

u e{ r ’ = r)= u e 


(4c, d) 


where Ri is the radius of the trough and W the width of the wheel. The Navier-Stokes Equations, continuity equation, 
and expression for the shear stress become: 5 


V 2 U g = 0 

dP _ pu 2 g 
dr r 

W J u g dr = constant 

gap 


(5a) 

(5b) 

(5c) 
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^ = n 


(5d) 


' du g Ug ' 

v dr r ) 


where P is the pressure in the gap. The solution to Eq. (7a), satisfying the boundary conditions (6a-d), is the 
following (subscript for k not shown for simplicity): 


Ug 


u^ a 


- Kj {kRj )/i ( kr ) + I x ( kR 2 )K X (kr) 
- K l ( kR 2 )/, (kR t ) + /, (kR 2 )K X ( kR l ) 


cos (kz) 


( 6 ) 


in which k assumes the values 


nn 


k ^ =k n) = ~^ (« = 1,3,5,...) 

W 


(7) 


to satisfy the boundary conditions (4c,d). In Eq. (6) 7) and K\ are modified Bessel functions of the first and second 
kinds. 

Now, it will always be true that kr » 1 for all r in the range R\ < r < R 2 , in which case 7) and K x approach their 
asymptotic limiting expressions: 


I\( kr ) 

Ki (kr) — > 


exp(A'r) 

yjlTtkr 

(8a) 

exp(-Ar) 

(8b) 

V 2nkr 


both of which are also valid for r = R\ and R 2 . Upon substituting (8a, b) into (6), one can rewrite the latter as 


Ug 


u y a pT sinh[/,(/e 2 -/•)] 

„ \ r sinh[A(i? 2 -i?,)] 


cos(A;z) 


Then, substituting Eq. (9) into (5d) yields for the shear stress anywhere in the gap 


(9) 


r 


re 



kg x cosh[A:(7? 2 


r)] 


3 

2 


\ r J 


sinh[A(i? 2 


sinh[A(i? 2 - R ] )] 


r)] 

— cos(Az) » 


(10) 


where T 0 is the stress in an equivalent flat-plate gap 
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flU 

gl 


(11a) 


r 


and 


Si — R 2 —R 1 


(lib) 


is the gap prior to the insertion of the sensor to be calibrated. (The negative sign in Eq. (10) is related to the 
definition of the positive radial direction.) At the boundary r = R h the condition u e = U is realized if the Fourier 
coefficients in Eq. (9) are chosen as 


2 


= (-D 


,,-if 2 y 


\ nn J 


(n = 1, 3, 5, ...) 


(12a) 


(12b) 


At the surface of the mounting plate ( r = R : ) the shear stress becomes 


T 


rd 



1 + y ( 2_V_ A# kg 1 cos (kz) 

2 „ l sinh(^g 1 ) n 


(n = 1,3,5,...) 


(13) 


The mean value of the shear stress across the width of the wheel becomes 


T rff 



jgi 1 

sinh(^gj) it 


{n= 1,3,5,...) 


(14) 


The factors in Eq. (14) represent the flat parallel-plate shear stress, the correction due to the curvature of the 
wheel, and the correction for the finite width of the wheel, respectively. Upon insertion of the test sensor into the 
trough an additional correction is needed to account for the flatness of the sensor. Then the expression for the shear 
stress can be written 


r,.,=-r 0 (l + C 1 )(l + C 2 )(l + C 3 ) 


( 15 ) 


where C 1 , C 2 ,and C 3 are the three respective corrections. 

B. Correction for the Curvature of the Wheel 

The correction follows from Eq. (13): 
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( 16 ) 


c, = E-.-IIl 

■y i? 2 2 /?, 

with the aid of Eq. (1 lb). Using values of the example at the end of Section II, the correction becomes 

1 1 £ 1 A - 6 

C, = — x— = -7.6 xl(TV- 0.0076%. 

1 2 0.5 


C. Correction for the Finite Width of the Wheel 

The correction is attributable to the terms containing kg, in Eq. (13). Without these terms the summation 
would have a value of A. Therefore 



' sinh(£gj ) n 


(17) 


The magnitude of the correction is shown as 
a function of g\!W in Fig. 4. If the wheel has a 
width of 0.0508 m and the gap remains 76 pm, 
then g\/W = 0.0015 and the correction amounts 
to C 2 = -6.6 x 10’ 4 or -0.066%. To achieve a 
correction not exceeding 0.01% would require a 
ratio g\/W = 0.0001, corresponding to an 
impractical wheel width of 0.76 m. 

D. Correction for Flatness of Test Sensor 

The sensor is mounted into a hole in the 
bottom surface of the trough with its rim flush 
at the intersection. Then the surface of the 
sensor will protrude into and thus narrow the 
gap, as illustrated in the front and top views 
shown in Fig. 5. However, there will be no 
surface discontinuities. The narrowing of the 
gap over the area of the sensor will lead to an 
increase in shear stress above that indicated in 
Eq. (11a). The correction will be based upon 
the mean value of the gap g(x) over the surface 
of the sensor. 

Let go, g i, and g(x) be the gap thicknesses at the center of the sensor, at the rim of the sensor, and at any 
distance x from the center of the sensor, respectively. Geometrical considerations reveal that 



Figure 4. Magnitude of the correction C 2 due to the 
finite width of the wheel. 


1^1 + g\)~~ a2 - [X +g(*)] “ _ * 2 - 1^1 + go ) 2 ( 18a ) 
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Figure 5. Detailed geometry of the sensor mount showing 
both front and top views. 


from which the variable gap g(x) is found: 

g(x) = (i ?2 ~ a' + X 2 ) 1 ' 2 -i?i (18b) 


The elemental area dA about a chord at the distance x from the center is 

dA = 2[a 2 -x 2 ) 1/2 d.v (19) 

Then the mean gap becomes 
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1 


(20) 


g = 


-(2 [fe-oW 

na J 



x 2 ) l,2 dx 


Under the assumption that a/R 2 «1, the integral can be evaluated in closed form, with the result: 


g = 



1 a 2 - x 2 

2 R, 


( a 


x 2 ) 1/2 dx = gl fl 

V 


3 a 2 ' 
8 g\ R l; 


(21) 


Upon inserting Eq. (21) into (11a), one finds the correction 


_ 3 a 2 

^3 


8 gl R 2 


(22) 


Using the parameter values of the preceding examples, together with a sensor having a radius of 0.00127 m (0.05 


„ J 

in.), one finds that the correction amounts to C, = — x 


0.00127 2 


= 0.016 or 1.6%. This correction, in 


8 76x10 6 x 0.5 

contrast to the others, leads to an increase in the shear stress. 

The small constriction in the flow channel caused by the presence of the sensor causes an acceleration of the 
flow into the gap. If the flow remains incompressible, then the invariance of the volume velocity causes Eq. (5c) to 
take the following form: 


-Ug 1 


u 0 = 


2g 


(23) 


the effect of which is accounted for in the correction C 3 . Unfortunately C 3 also shows a limiting factor for this 
design. As the radius of the sensor increases the correction increases by the square. For the design example if we 
limit C 3 to 0.1 that imposes a limit on the sensor radius of 0.0032m (0.13 in.). If a larger gap size can be used to 
generate the needed shear stress then the size of the sensor may increase and maintain the value of C 3 . 

E. Stability of the Flow 

The rotary flow channel represents a special case of Couette flow, which displays a multitude of instability 
modes even at low Reynolds numbers. The criterion for the onset of instability is based on the Taylor number 6 : 


Ta = 


2p 2 orR x gl 

. 2 




(24) 


A Taylor number not exceeding the critical value of 1708 is said to insure a stable flow. Upon inserting above 
parameter values and those found in Table I, one finds the Taylor number (in neon): 
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= 181 


Ta = 


2 x 0.839 2 x 754 2 x 0.5 x (76 x 10~ 6 ) 3 
(3.1 lx 1CT 5 ) 2 


Since the Taylor number is far below critical, the rotary channel flow is expected to be stable. 

F. Centrifugal Pressure on the Sensor 

The centrifugal pressure on the sensor follows directly from Eq. (5b): 



u 2 g dr 

r 


Upon inserting the value of u g in the limit of small kg t 


U{R 2 -r ) 

Si 


(25) 


(26) 


one finds 


P(R 2 ) = -p 


r u \ 

2 


( rA 


Rl In 

2 

U.j 


L 



+ f(*> 


-3 R 2 ) 


(27) 


With the parameter values used above, the centrifugal pressure exerted on the sensor amounts to 2.35 atm. This will 
not load the sensor if the sensing element is back vented. Further, the Knudsen number, defined here as the mean 
free path over the thickness of the gap, is 9.4 x 1 0“ 8 / 76 x 10' 6 = 0.0012, indicating that the flow is well within the 
domain of continuum theory. 


IV, IV. Design Review and Possible Extensions 

By choosing a single value for the radius of the wheel, R u we obtain a simple correlation between the intended 
gap (g!), peripheral speed of the wheel (U), and the generated shear stress as shown in Eq. 3. Of more interest is a 
design guideline for range of shear stress in air and neon for different gap sizes. Table 2 and Table 3 show this 


Table 2: Achievable shear stress versus gap size in 
air at 20°C for /?i=0.5m and minimum ro=60 rpm, 
the mach number limited shear stress at m= 1899.74 
rpm, and maximum shear stress at r«=5500 rpm. 


Gap, pm 

^ mini P& 

^machi P^ 

Trnax, Pa 

400 

0.14 

4.55 

13.18 

300 

0.19 

6.07 

17.57 

200 

0.29 

9.10 

26.35 

150 

0.38 

12.14 

35.13 

100 

0.57 

18.20 

52.70 

75 

0.77 

24.27 

70.27 

50 

1.15 

36.41 

105.4 


Table 3: Achievable shear stress versus gap size in 
neon at 20°C for I?i=0.5m and minimum co = 60 rpm, 
the mach number limited shear stress at ro=2486.83 
rpm, and maximum shear stress at ro=7200 rpm. 


Gap, pm 

t mini P& 

^ machi 

^ maxi P^ 

400 

0.24 

10.12 

29.31 

300 

0.33 

13.50 

39.08 

200 

0.49 

20.25 

58.62 

150 

0.65 

27.00 

78.16 

100 

0.98 

40.50 

117.24 

75 

1.30 

53.99 

156.33 

50 

1.95 

80.99 

234.49 
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information for air and neon respectively. The minimum shear stress is based on co = 60 rpm and scales linearly for 
decreased values of to. The shear stress listed for z mac h corresponds with a linear velocity limited to M= 0.29. The 
maximum shear stress corresponds to a linear velocity where M= 0.84. 

These tables are useful for quickly determining the gap size needed to generate the shear stress range of interest. 
Comparing these values with the shear stress obtainable using a flow channel shows that the flow channel generates 
a higher shear stress for a given height. One advantage with the rotary flow channel is while it needs to be 
completely enclosed, there is no exhaust to the ambient environment required. This simplifies the use of neon for the 
higher viscosity to obtain higher shear stress. 

The rotary flow channel, as well as the tape based design; both have the potential to generate a small dynamic 
shear stress about a mean value. This dynamic component while small compared to the maximum shear stress 
attainable would be achievable with a wide range of mean values at a range of frequencies dependent on the details 
of implementation. The method to do this would be to very precisely move or vibrate the trough at a known 
magnitude and frequency. The analysis of this flow would be more complex than the analysis shown here since the 
gap would not be constant except at the initial position of the trough as analyzed previously. 

V. Conclusions 

The rotary flow channel has many desirable characteristics for generating a known shear stress with the potential 
for providing a dynamic shear stress capability as well. The corrections needed for a practical sized device are small 
and easily determined. To keep the corrections small does impose a very stringent requirement on the sensor size but 
for MEMS devices this should not be insurmountable and may even be advantageous in that a flush mounted device 
may be able to be pre-installed in a plug type test article that is then mounted with only a portion of the sensor die 
exposed to the calibration flow. Additional analysis needs to be done with particular emphasis on wall slip 
conditions. Since this is not a pressure driven flow compressibility should not be an issue with the possible 
exception of the region at the sensor installation due to the non-uniformity in the gap in this region. 
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